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On  the  Geometry  of  Sculptured 
Surface  Machining 


Johannes  Wallner  and  Helmut  Pottmann 


Abstract.  We  present  geometric  aspects  of  sculptured  surface  machin- 
ing. Several  possible  configuration  manifolds  of  tool  positions  relative  to  a 
workpiece  are  investigated  under  different  aspects:  the  degree  of  freedom 
of  the  motion  of  the  tool,  the  correspondence  between  the  contact  point 
and  the  tool  position,  and  the  presence  or  absence  of  unwanted  collisions 
between  tool  and  workpiece. 


§1.  Introduction 

In  the  past  decades,  strong  research  efforts  have  been  devoted  to  developing 
the  mathematical  fundamentals  and  efficient  algorithms  for  the  representation 
of  free-form  surfaces  in  CAD/CAM  systems.  However,  just  a few  contributions 
address  manufacturing  of  sculptured  surfaces,  although  there  are  appealing 
and  practically  important  open  problems  in  this  area. 

Geometrical  problems  in  this  area  include  the  following:  If  two  surfaces 
touch  each  other  at  a point,  such  as  a milling-tool  and  a free-form  surface 
which  is  to  be  manufactured,  does  the  curvature  of  the  surfaces  force  them 
to  intersect  arbitrarily  near  the  contact  point  (=  the  local  collision  problem)? 
Given  a free-form  surface  and  a milling  tool,  is  the  tool  able  to  move  such 
that  its  envelope  during  the  motion  is  the  free-form  surface?  Is  the  tool  able 
to  do  this  while  moving  only  by  translations  (3-axis  milling),  or  do  we  need 
more  flexibility  (5-axis  milling)?  Which  relative  tool  position  achieves  best 
surface  quality  (=tool  positioning)?  How  can  we  decompose  the  theoretically 
two-parameter  motion  of  the  milling  tool  by  a series  of  one-parameter  motions 
such  that  e.g.  manufacturing  time  is  minimal  (=motion  planning)? 

A survey  of  mathematical  fundamentals  on  NC  machining  of  sculptured 
surfaces  is  given  in  [1,14,19].  Tool  selection,  motion  planning  and  local  inter- 
ference checking  for  3-axis  and  5-axis  machining  has  been  studied  in  [2,9,10,11, 
12,13,25].  Three-axis  machining  (general  offsets  and  Minkowski  addition)  are 
considered  by  [15,21,23].  A configuration  space  has  been  defined  in  [3,24], 
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There  is  a variety  of  contributions  using  Computer  Graphics  techniques 
such  as  visibility  algorithms  [5,7,22].  There  the  cutter  shaft  is  shrunk  to  its 
axis,  and  simultaneously  the  design  surface  is  replaced  by  an  appropriate 
offset  surface.  Collision  checking  is  thus  transformed  to  a visibility  test  (5- 
axis  machining  requires  a modification  of  the  method). 

This  paper  is  organized  as  follows:  In  Section  2 we  briefly  show  some 
concepts  of  elementary  differential  geometry.  Section  3 sums  up  some  results 
concerning  the  local  contact  situation.  Section  4 investigates  configuration 
manifolds  of  motions  constrained  in  several  ways  and  describes  the  possible 
infinitesimal  motions  of  a milling-tool  in  a contact  position.  In  Section  5 we 
study  the  dependency  of  the  contact  point  of  the  contact  position.  Finally 
Section  6 features  global  statements  about  the  absence  of  unwanted  collisions 
under  certain  circumstances. 

§2.  Differential  Geometry 

We  first  give  a short  description  of  some  aspects  of  curvature  theory  of  2- 
surfaces  in  Euclidean  three-space  (cf.  [4,17,20]). 

2.1.  Oriented  surfaces  and  their  first  and  second  fundamental  forms 

Consider  a regular  smooth  surface  given  by  the  parametrization  / = f(u), 
where  x = (xi , a?2»  *3)  is  a point  in  Euclidean  1R3  and  u = (1*1,112)  ranges 
in  some  open  planar  domain  D.  We  assume  that  / is  twice  continuously 
differentiable  in  order  to  be  able  to  define  curvatures. 

The  differential  df  of  / maps  a tangent  vector  v = (iq,  V2)  attached  to  the 
point  u,  to  the  vector  duf(v)  — -^\t=of{u+tv),  which  is  computed  by  duf(v)  = 
xti(u)v\  + x, 2(11)1)2,  where  the  symbols  x,i  and  x,2  mean  differentiation  with 
respect  to  the  first  and  second  variable. 

The  function  n = (x:i  x xi2)/||xii  x x^H  is  the  surface  unit  normal  vector. 
The  symmetric  bilinear  forms  gu(v,w)  = duf(v)  ■ duf(w)  and  hu(v,w)  = 
—dn(v)  ■ w are  called  the  first  and  second  fundamental  forms  of  /.  With  g,j  = 
Xj-Xj  and  hjj  = —nyxtj  — n-xy  we  have  g(v,  w)  = and  h(v,w)  — 

y ) hij  Vi  Wj . 

2.2.  The  Dupin  indicatrix  and  Meusnier’s  theorem 

If  / is  the  parametrization  of  a surface,  and  u = u(t)  is  a curve  in  its  parameter 
domain  D , then  e(<)  = f(u(t))  is  a curve  contained  in  the  surface  f(D).  The 
curve’s  tangent  c(t)  = df(u(t))  is  contained  in  the  surface’s  tangent  plane.  Its 
second  derivative  vector  is  split  into  three  components: 

c = ||c||  («„n  + qc  + Kgb) , 

where  n is  the  normal  vector  evaluated  at  u(t)  and  b is  the  curve’s  normal  in 
the  tangent  plane.  The  coefficients  Kg  and  Kn  are  the  geodesic  curvature  and 
normal  curvature  of  the  curve,  respectively.  The  following  theorem  states  the 
perhaps  unexpected  fact  that  the  normal  curvature  is  dependent  only  on  the 
direction  c,  and  we  can  therefore  speak  of  the  normal  curvature  of  a surface 
tangent. 
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Theorem  1.  (Meusnier)  The  normal  curvature  of  the  curve  c is  computed 
by  Kn  = nn(u)  = h(u,u)/g(u,u). 

Evaluate  the  matrix  product  (<?ij)_1hjj  at  and  compute  its  two  lin- 

early independent  eigenvectors  v',v"  which  correspond  to  eigenvalues  Hi,  k2. 
Then  df(v'),  df(v")  are  orthogonal  and  define  the  two  principal  surface  tangents 
at  u. 

Theorem  2.  (Euler)  Assume  that  w',w"  are  unit  vectors  parallel  to  the 
principal  surface  tangents  at  u.  If  df(v)  = cos  <f>w'  + sin  <f>w",  then  Kn(v)  = 
cos2  <f>  ■ Ki  + sin2  <j)  ■ K2-  The  polar  diagrams  of  1 / \JiTf.  and  l/y/—nn  in  the 
tangent  plane  (the  oriented  Dupin  indicatrices  i+,  i-)  are  possibly  void  or 
singular  conic  sections  centered  in  the  origin. 

Surface  points  are  called  elliptic,  if  ki«2  > 0,  hyperbolic  if  ki k2  < 0,  flat  if 
ki  = K2  = 0,  and  parabolic  in  the  remaining  cases. 

2.3.  Euclidean  displacements  and  infinitesimal  motions 

A Euclidean  displacement  g : x € R3  h->  g(x)  € 1R3  may  be  written  in  the 
form  x M ■ x + v,  where  M is  an  orthogonal  matrix  of  determinant  1.  A 
one-parameter  family  g(t)  = (v(t),M(t))  of  Euclidean  displacements,  (=  a 
path  of  Euclidean  motions,  or  a smooth  Euclidean  motion)  has  in  all  of  its 
instants  an  infinitesimal  motion,  which  is  determined  by  the  velocity  vectors 
d(g(t))(x)/dt  = v(t)  + M(t)-x  of  all  points.  If  an  infinitesimal  motion  coincides 
with  the  velocity  field  of  a smooth  rotation  about  an  axis,  it  is  called  an 
infinitesimal  rotation.  The  definition  of  infinitesimal  translation  and  infinitesimal 
helical  motion  is  analogous. 

There  is  a linear  space  of  infinitesimal  motions.  It  is  further  well  known 
that  all  infinitesimal  motions  (i.e,  all  velocity  fields  of  smooth  motions)  can 
be  written  in  the  form 

x = c + c x x. 

The  condition  c-c  = 0 characterizes  infinitesimal  rotations,  c = 0 characterizes 
infinitesimal  translations,  and  c ^ 0,  c-c  ^ 0 characterizes  infinitesimal  helical 
motions.  We  briefly  write  (c,  c)  6 1R6  to  denote  an  infinitesimal  motion.  An 
infinitesimal  rotation  whose  axis  is  a + [6]  then  has  the  form  A (b,  a x b)  (the 
symbol  [6]  denotes  all  multiples  of  the  vector  b ).  The  translation  x = c has  the 
form  (0,  c).  If  we  have  to  consider  the  coordinates  of  c,  c in  some  coordinate 
system,  we  always  write  c = (c0i, c02, C03)  and  c = (c23, C31, ci2). 

2.4.  Ruled  surfaces 

If  p(t),  v(t)  are  two  curves  with  v / 0,  then  f(u i,w2)  = p(«i)  + it2v(ui)  is 
a surface  whose  parameter  lines  u\  — const  are  straight  lines,  which  will  be 
denoted  by  l(u\).  Such  a surface  is  called  a ruled  surface. 

We  need  the  following  well  known  results  concerning  the  first  order  differ- 
ential properties  of  ruled  surfaces  in  Euclidean  space:  there  is  an  orthonormal 
frame  ( q ; ei,  e2,  e3),  dependent  on  a parameter  t,  and  a smooth  function  Ui (t) 
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such  that  q(t)  is  on  the  ruling  l(ui(t)),  e\ (t)  is  parallel  to  it,  and  e-i (t)  is 
tangent  to  the  surface  at  q(t ),  and  which  can  be  chosen  such  that  its  infinites- 
imal motion  at  t = 0 is  a helical  motion  (c,  c)  which  can  be  computed  by  the 
following  formulae: 


p p x v s := 


lb  x Pll 


2 (det (p,p,p)p  + det(p,p,p)p  + (p  ■ p)(p  x p)) 


b = 


(p  x v)  x v 


II (p  xv)  x n|| 


( pp)p 2 


(=0 


(p  x p)2 


c = <56,  c—b—cx  s(0). 


t=0 


The  first  order  differential  invariant  6 is  called  the  distribution  parameter  of 
the  ruled  surface,  and  the  point  s(t)  is  called  its  striction  point. 


2.5.  Line  congruences 

A smooth  line  congruence  /C  is  a smooth  two-parameter  family  of  straight 
lines  l(u\,U2)  in  Euclidean  three-space.  It  may  be  parametrized  by  two  ‘sur- 
faces’ p(ui,U2)  and  v(ux,U2),  where  v(u\,U2)  / 0.  The  line  l(u  1,1/2)  then  is 
p(u\ , 112)  + b(ui,«2)].  The  choice  of  a curve  (ui(t),U2(<))  in  the  parameter 
domain  gives  a ruled  surface  f(ux(t), u2(t)). 

Definition.  A line  congruence  fC  is  regular  at  a line  l if  the  six-tuples  (p,pxv), 
(p,iiP  X t),i  + p,  1 x v),  (p,2,P  x «,2  + p, 2 X v)  are  linearly  independent. 

The  meaning  of  this  definition  is  that  the  lines  of  the  congruence  actually 
change  infinitesimally  if  we  move  infinitesimally  in  the  parameter  domain. 
We  will  always  assume  that  1C  is  regular. 

It  is  well  known  that  the  infinitesimal  properties  of  first  order  of  a line 
l within  K.  are  like  those  of  a linear  line  congruence  1C',  which  is  called  the 
tangent  to  1C.  There  are  the  following  possibilities  for  K'\ 

• K'  is  the  set  of  lines  which  intersect  two  lines  k' , k"  (=  a hyperbolic  linear 
congruence  with  axes  k',  k"). 

• 1C'  is  the  set  of  lines  whose  complex  extensions  intersect  two  conjugate 
complex  lines  k' , k"  (=  an  elliptic  linear  congruence).  This  set  of  lines 
is  the  affine  image  of  the  set  of  lines  which  join  the  points  (x,  y,  0)  and 
(a:  cos  <j)  — y sin  <f>,  x sin  cf>  + y cos  <j>,  1),  0 < <j)  < 27r. 

• K!  is  the  set  of  lines  tangent  to  a ruled  quadric  in  the  points  of  one  of  its 
rulings  k.  (=  a parabolic  linear  congruence).  The  line  k is  also  called  the 
axis  of  1C’. 

• 1C'  is  a bundle  of  lines. 

In  all  cases,  lines  at  infinity  are  allowed  (but  in  the  elliptic  case  they  do  not 
occur).  These  four  types  of  lines  l in  a congruence  are  accordingly  called 
hyperbolic,  elliptic,  parabolic,  and  degenerate. 

In  the  hyperbolic  and  parabolic  case,  the  points  of  l contained  in  an  axis 
are  called  focal  points,  in  the  degerenate  case  the  bundle  vertex  is  also  called 
a focal  point. 
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There  is  also  a connection  between  infinitesimal  motions  and  linear  line 
congruences:  a four-dimensional  subspace  of  the  space  of  infinitesimal  mo- 
tions always  contains  infinitesimal  rotations  or  translations,  whose  axes  form 
a linear  line  congruence  (here  we  assign  to  the  infinitesimal  translation  x = c 
the  axis  at  infinity  which  is  orthogonal  to  c). 

2.5  Submanifolds 

As  the  concepts  listed  here  are  essential  only  for  the  proofs  of  our  results, 
we  give  only  a brief  summary:  We  assume  that  the  reader  is  familiar  with 
the  concept  of  a smooth  n-dimensional  manifold  N (an  n- manifold).  An  em- 
bedded smooth  m-submanifold  M of  N is  characterized  by  the  existence,  for 
p e M,  of  a local  diffeomorphism  which  transforms  an  iV-neighbourhood  U of 
p to  ]Rn  and  U fl  M to  Htm  C IR".  An  immersion  is  a smooth  mapping  whose 
differential  is  one-to-one  (but  not  necessarily  onto).  Then  locally  the  immer- 
sion is  also  one-to-one.  An  immersed  /c-submanifold  is  the  image  of  a smooth 
fc-manifold  under  an  immersion.  The  difference  between  embedded  and  im- 
mersed submanifolds  is  therefore  that  the  latter  may  have  ‘self-intersections’, 
but  neither  are  allowed  to  have  ‘singularities’. 

An  embedded  submanifold  Mi  and  an  immersed  submanifold  M2  of  N 
are  transverse  (we  write  Mi  ifl  M2),  if  for  all  points  p € Mx  n M2  the  tangent 
spaces  TpMi,  TVM2  span  TPN.  Then  Mi  fl  M2  is  an  immersed  (dim Mi  + 
dim  M2  — dim  iV)-dimensional  submanifold  of  N,  whose  tangent  space  equals 
TPM\  fl  TpM-i- 

If  a smooth  mapping  <f  of  a smooth  m-manifold  M into  a smooth  71- 
manifold  N has  constant  rank  r (i.e.,  at  all  points  its  differential’s  rank  as  of 
a linear  mapping  equals  r),  then  0(M)  is  a smooth  immersed  r-dimensional 
submanifold  of  N. 


§3.  Local  Contact  Situation 

If  a body  is  bounded  by  a smooth  surface  /,  this  surface  has  an  inside  and 
an  outside.  The  unit  surface  normals  can  point  to  either  side,  depending  on 
the  parametrization.  If  two  bodies,  which  are  bounded  by  smooth  surfaces 
/",  touch  each  other,  the  curvatures  of  /',  /"  give  information  whether  they 
intersect  locally  or  not. 

It  should  be  remarked  that  some  methods  proposed  in  the  literature  (cf. 
[10,12])  for  avoiding  local  intersections  are  only  approximations,  and  one  can 
find  surfaces  where  they  won’t  work.  Also  it  is  important  to  note  that  the 
presence  or  absence  of  local  intersections  is  completely  independent  of  the 
actual  motion  of  the  two  bodies. 

Definition.  The  interior  int(i)  of  a conic  i centered  in  the  origin  is  void  if  i 
is  void,  and  otherwise  is  the  connected  component  of  IR2  \ i which  contains  0 
and  whose  boundary  is  i.  Its  exterior  ext  (?)  is  the  complement  ofint(i)  U i. 

We  assume  that  /',  f"  are  parametrized  such  that  their  unit  normal  vectors 
in  the  common  point  coincide,  and  that  this  common  normal  vector  points  to 
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the  inside  for  /'  and  to  the  outside  for  /".  We  further  assume  that  not  all 
four  principal  curvatures  are  zero.  We  consider  the  Dupin  indicatrices  i'+ , i'_ , 
i'_[,  and  i"_  for  /'  and  f",  respectively. 

Theorem  3.  Under  these  assumptions,  the  two  bodies  in  question  intersect 
locally,  if  one  of  the  intersections  z"  D int(z'+),  i'L  fl  ext(?'_)  is  not  void.  They 
do  not  intersect,  if  i"  fl  i'+,  i'\.  fl  int(z'+),  i"_  fi  i'_,  i'L  fl  ext (i'_)  are  void. 

Proof:  (Sketch)  The  theorem  follows  from  the  fact  that  a twice  continuously 
differentiable  surface  may  be  approximated  of  second  order  by  the  graph  of 
a quadratic  function,  whose  contour  lines  are  scaled  versions  of  the  Dupin 
indicatrices.  □ 

Note  that  the  theorem  says  nothing  about  the  cases  that  all  principal 
curvatures  are  zero,  the  indicatrices  touch  each  other  in  two  opposite  points, 
or  even  coincide  (cf.  Fig.  1,  right).  In  that  case,  second  derivatives  are  not 
sufficient  to  decide  if  there  are  local  intersections.  In  practice,  this  does  not 
matter  very  much  because  the  only  case  that  is  likely  to  occur  with  nonzero 
probability  is  that  of  a flat  end  mill  shaping  a planar  surface,  which  does  not 
have  self-intersections. 

In  [8,16,23]  it  is  shown  how  to  define  indicatrices  in  the  case  of  piecewise 
curvature-continuous  surfaces.  The  theorem  is  valid  also  in  this  more  general 
case. 


§4.  Configuration  Manifolds  and  their  Tangent  Spaces 

For  many  problems  concerning  the  milling  of  free-form  surfaces,  it  is  important 
to  know  the  degree  of  freedom  of  a motion  constrained  in  various  ways.  Typical 
constraints  are:  Motion  by  translations  such  that  a surface  remains  in  contact 
with  another  surface  (3-axis  milling),  motion  such  that  a surface  remains  in 
contact  with  another  surface  (possible  set  of  tool  positions  in  5-axis  milling), 
motion  such  that  a milling  tool  remains  in  contact  with  a surface  and  its  axis 
is  contained  in  some  previously  prescribed  line  congruence  (a  possible  way  to 
do  5-axis  milling).  These  topics  will  be  discussed  in  Subsections  4. 1-4.3. 

4.1  Translational  motions  constrained  by  surface-surface  contact 

Consider  two  surfaces  /"  which  have  a common  point  p = f'(u')  — f"(u ") 
and  share  a common  unit  surface  normal  n'  = n"  there.  Imagine  the  first 
surface  moving  by  translations  such  that  it  always  touches  the  second  surface 
(i.e.,  at  every  instant  t there  are  a translation  vector  a(t)  and  parameter  values 
u'(t)  and  u"(t)  such  that  f'(u'(t))+a(t)  = /"(«"(<))).  One  would  expect  that 
this  motion  has  two  degrees  of  freedom,  if  we  do  not  count  intersections  of  the 
surfaces. 

If  q(xi,xi)  = YL  ajjXiXj  is  a bivariate  homogeneous  quadratic  polynomial 
in  the  variables  aq,  aq,  we  call  the  rank  of  the  (2  x 2)-matrix  atJ  the  rank  of  q. 
The  zero  set  of  q consists  of  the  entire  plane  in  the  case  of  zero  rank,  of  one 
line  if  the  rank  is  one,  and  of  two  real  or  two  conjugate  complex  lines  if  the 
rank  is  two. 
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q’-q”=0 

q’-q"=0 

Fig.  1.  Indicatrices  q — 1,  q"  = 1 and  zero  set  of  q — q" . Left:  rk (q'  — q")  = 2, 
Right:  rk (q1  — q)  — 1.  . 

We  have  the  following  theorem  (cf.  [8,16,23]),  which  is  valid  for  all  contact 
situations  where  at  least  one  surface  has  no  parabolic  or  flat  point,  or  both 
have  parabolic  points  but  the  principal  tangents  do  not  coincide. 

Theorem  4.  Write  the  equation  of  the  oriented  Dupin  indicatrices  of  f , /" 
at  the  contact  point  in  the  form  q'(x i,  *2)  = 1 and  q"{xi,x 2)  = 1,  where  x\ , x? 
are  Cartesian  coordinates  in  the  tangent  plane,  and  q' , q"  are  bivariate  homo- 
geneous quadratic  polynomials  in  X\,X2-  If  the  condition  stated  immediately 
before  this  theorem  is  satisfied,  then  the  rank  of  q'  — q"  gives  the  infinitesimal 
degree  of  freedom  of  translational  motions  constrained  by  the  contact  of  f, 
f"  (see  Fig.  1). 

Proposition.  If,  under  the  assumptions  of  Th.  4,  the  infinitesimal  degree  of 
freedom  is  two,  then  so  is  the  local  degree  of  freedom. 


4.2  Motions  constrained  by  surface-surface  contact 

Definition.  The  set  of  proper  Euclidean  motions  which  transforms  a surface 
f such  that  it  touches  a surface  f"  is  called  the  configuration  space  C = 
C(f,  f")  of  surface-surface  contact. 

Clearly,  a position  g € C is  not  determined  by  the  contact  points  alone, 
because  we  still  may  rotate  /'  about  the  contact  normal.  But  if  we  prescribe 
a unit  tangent  vector  ( p'-,w ')  of  /'  and  (p';  w")  of  f"  (which  means  p'  — 
f'(ui,u 2)  and  w'  = df'(v')  with  ||w'||  = 1,  and  the  same  for  (p” ,u>")),  then 
there  is  a unique  Euclidean  motion  g € C which  maps  not  only  the  point  p' 
onto  p",  but  also  the  tangent  vector  w'  onto  w".  If  we  rotate  both  w' , w" 
about  an  angle  (j),  this  leads  to  the  same  p € C,  so  (7  is  a smooth  image  of  the 
factor  manifold  C of  such  equivalence  classes  of  unit  tangent  vectors,  which  is 
five-dimensional. 

This  shows  that  we  may  expect  five  degrees  of  freedom,  if  /'  moves  under 
the  constraint  that  it  touches  f"  in  some  point.  The  following  theorem  is 
given  in  [18,24]: 
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Theorem  5.  In  the  notation  of  Th.  4,  the  motion  of  a surface  under  the 
single  constraint  that  it  touches  a second  one,  has  five  degrees  of  freedom  if 
q'  — q"  has  rank  two.  If  that  is  the  case  for  all  possible  contact  points,  then 
the  configuration  space  is  an  immersed  five-dimensional  submanifold  of  the 
motion  group. 

Note  that  the  theorem  is  valid  without  the  additional  assumption  made  in 
Th.  4 in  case  some  principal  curvatures  are  zero.  It  is  easy  to  describe  the 
tangent  space  TgC  of  the  configuration  manifold  at  a contact  position  g £ 
C C G ( to  be  precise,  the  tangent  space  of  the  immersion  described  at  the 
beginning  of  Sect.  4.2).  The  proof  of  the  following  proposition  can  be  found 
in  [18]. 

Proposition.  The  linear  space  TgC  of  infinitesimal  motions  which  belong  to 
paths  in  C is  five-dimensional,  and  the  axes  of  its  infinitesimal  rotations  are 
the  lines  which  intersect  the  contact  normal. 

In  a Cartesian  coordinate  system  whose  origin  is  the  contact  point  and  whose 
a:3-axis  is  the  contact  normal,  TgC  has  the  equation  C03  = 0. 

4.3.  Motions  constrained  by  congruences 

Here  we  consider  the  motion  of  a rigid  body  E such  that  a line  a of  E is 
contained  in  a congruence  K.  This  subset  of  the  group  G of  proper  Eu- 
clidean motions  will  be  denoted  by  K.  We  assume  that  K.  is  parametrized  by 
l(ui,u2)  = p{u  i,u2)  + [v(ui,u2)]. 

If  g € K,  then  it  is  obvious  that  both  t o g and  po  g are  in  K,  if  r is  a 
translation  parallel  to  g{a ) and  p is  a rotation  with  axis  g(a). 

Lemma.  Assume  that  1C  is  a smooth  line  congruence  which  is  regular  at  g(a). 
Then  K is  a four-dimensional  smooth  submanifold  of  G in  a neighbourhood 
ofg. 

Proof:  (Sketch)  Let  g € K be  a position  of  E such  that  g(a)  e fC.  We 
compute  K’s  tangent  space  TgK  of  infinitesimal  motions  at  g:  Consider  a 
curve  (ui(t),v.2(t))  in  1C’ s parameter  domain  such  that  l(t)  = l(ui(t),U2(t)) 
is  a ruled  surface  within  K with  1(0)  = g(a).  The  helical  motion  described 
in  Sect.  2.4  is  tangent  to  K.  If  we  choose  two  such  curves  with  linearly 
independent  tangent  vectors,  this  gives  two  linearly  independent  infinitesimal 
motions  of  TgK,  if  K is  regular  at  g(a)  (the  proof  of  this  is  left  as  an  exercise 
to  the  reader). 

Obviously  all  infinitesimal  translations  parallel  to  g(a)  are  in  TgK,  and 
so  are  the  infinitesimal  rotations  with  axis  g(a).  If  g(a)  is  the  line  p+  [v],  then 
the  former  is  described  by  the  six-tuple  (0,  v)  and  the  latter  by  (v,p  x v ).  Now 
TgK  is  the  linear  span  of  these  four  infinitesimal  motions.  □ 

In  case  that  g(a)  is  a hyperbolic  line  of  1C,  there  is  a simple  geometric 
characterization  of  the  infinitesimal  rotations  of  TgK: 
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Fig.  2.  Axes  of  infinitesimal  rotations  in  K. 

Proposition.  If  1C  is  a hyperbolic  linear  congruence  tangent  to  1C  at  the  line 
g(a)  and  has  axes  k',  k" , then  the  axes  of  infinitesimal  rotations  in  TgK  form 
a hyperbolic  linear  congruence  C,  whose  axes  b',  b"  are  incident  with  g(a)C\k", 
g(a)  n k',  respecively,  are  orthogonal  to  g(a),  and  are  such  that  k',b',g(a)  as 
well  as  k",b",g(a)  are  coplanar  (see  Fig.  2). 

Proof:  It  is  easy  to  see  that  the  lines  a',  a",  incident  with  k1  Ci  g(a),  k"  fi 
g(a)  and  orthogonal  to  k\  k" , respectively,  are  axes  of  infinitesimal  rotations 
contained  in  TgK.  The  axis  of  the  infinitesimal  translations  along  g(a)  (which 
is  the  line  at  infinity  orthogonal  to  g(a)),  g(a)  itself,  and  o',  a"  intersect  both 
b b" . We  already  know  that  the  set  of  axes  is  a linear  congruence;  because 
of  these  four  intersections  the  lines  b\  b"  are  necessarily  the  axes  of  C.  □ 

The  following  is  used  later: 

Proposition.  A point  of  g(a)  which  contains  two  different  axes  of  infinitesi- 
mal rotations  ofTgK  must  be  a focal  point  of  g{a). 

Proof:  At  least  one  of  the  two  axes  is  not  g(a)  itself,  and  the  rotation  about 
this  axis  transforms  g(a)  into  a line  which  intersects  g(a).  Looking  at  the  list 
of  linear  tangent  congruences  in  Sect.  2.4  shows  that  this  is  only  possible  in  a 
focal  point.  □ 

4.4.  Multiple  constraints 

Assume  that  a rigid  body  E,  bounded  by  a smooth  surface  /',  moves  such 
that  /'  remains  in  contact  with  a surface  and  that  in  addition  a line  a of 
E is  contained  in  a smooth  line  congruence  K.. 

With  C being  the  configuration  space  of  surface-surface  contact  (see 
Sect.  4.2)  and  AT  as  in  Sect.  4.3  the  set  of  possible  positions  g of  E is  given 
by  the  intersection  C fl  K.  The  following  theorem  shows  under  what  circum- 
stances Cf\K  is  actually  a smooth  three-dimensional  submanifold  of  positions, 
as  is  to  be  expected  when  comparing  dimensions: 

Theorem  6.  IfTgK  is  not  contained  in  TgC,  then  C Pi  K is  a three-dimensi- 
onal immersed  submanifold  in  a neighbourhood  of  g.  This  is  always  the  case 
if  g(a)  is  not  parallel  to  the  contact  tangent  plane. 

Proof:  We  have  to  show  that  C ifl  K.  Because  dimTgC  = 5,  this  is  always 
the  case  if  TgK  is  not  contained  in  TgC.  If  the  contact  tangent  plane  is  not 
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parallel  to  /,  then  any  small  translation  parallel  to  / leaves  C and  therefore  in 
this  case  C ffl  K,  and  dim((7  fl  K)  = dim  C + dim  K — dim  G = 3.  □ 

If  E has  rotational  symmetry,  and  the  line  g(a)  which  is  forced  to  belong  to 
K,  is  its  axis,  then  we  can  say  more  about  the  set  Cf\K  of  admissible  positions. 
Clearly,  any  rotation  about  the  axis  does  not  change  E,  so  g G C fl  K implies 
that  p o g g C fl  K whenever  p is  such  a rotation. 

We  choose  a reference  point  p on  the  axis,  and  look  at  the  three-parameter 
family  of  its  positions,  which  is  actually  only  a two-parameter  family: 

Proposition.  If  E has  rotational  symmetry,  a reference  point  p in  E ’s  axis  a 
traces  out  a regular  two-surface  while  undergoing  all  transformations  ofCfiK, 
provided  that  p is  never  a focal  point  of  g(a).  This  path  surface  is  transverse 
to  l if  l is  not  parallel  to  the  contact  tangent  plane. 

Proof:  We  consider  the  mapping  <j> : g i->  g(p)  of  Cfl  K to  1R3.  All  infinitesi- 
mal rotations  about  g(a)  assign  zero  velocity  to  p,  which  implies  that  the  rank 
of  (f>  is  not  greater  than  two.  If  p has  zero  velocity  also  for  other  infinitesimal 
motions  of  Tg(C  fl  K),  these  must  be  infinitesimal  rotations,  and  we  can  use 
the  proposition  at  the  end  of  Sect.  4.3  to  conclude  that  the  rank  of  tj>  is  indeed 
two,  and  its  image  a smooth  2-surface.  If  l and  the  contact  tangent  plane  are 
not  parallel,  then  no  infinitesimal  translation  parallel  to  / is  in  C Pi  K , and  the 
path  surface  cannot  be  parallel  to  l.  □ 

§5.  Movement  of  the  Contact  Point 

It  is  important  to  study  the  dependency  of  the  contact  point  on  the  contact 
positions.  This  will  be  done  for  two  different  types  of  constraints. 

5.1.  Motions  constrained  by  surface  contact 

In  Sect.  4.2  we  stated  that  the  configuration  space  C(/',  f")  of  surface-surface 
contact  is  an  immersed  image  of  C as  described  in  Sect.  4.2,  if  the  Dupin 
indicatrices  of  /',  f"  fulfill  a certain  condition  (Th.  5). 

If  this  is  the  case,  then  there  is  a local  inverse  C — > M,  and  so  the  contact 
point  depends  on  the  positions  g £ C in  a smooth  way  locally.  This  (local) 
mapping  will  be  denoted  by  i>. 

As  dim  C = 5 and  the  contact  point  varies  in  a 2-surface,  there  is  a 3- 
dimensional  kernel  subspace  ker  dgip  C TgC  of  infinitesimal  motions  which 
do  not  (infinitesimally)  change  the  contact  point.  The  following  is  an  easy 
exercise  in  differential  geometry: 

Lemma.  If  f,  f"  are  two  surfaces  having  contact  at  a point  p,  n is  the  unit 
normal  vector  of  f,  and  wi,  u>2  are  principal  tangent  vectors,  corresponding 
to  curvatures  k\,  ft 2,  then  all  infinitesimal  rotations  about  the  axes  p + [n], 
P + ti/ki  + [u)2j\  p + n//t2  + [wi]  are  contained  in  the  tangent  space  TgC  of  the 
configuration  manifold,  and  do  not  (infinitesimally)  change  the  contact  point 
on  f". 

If  a principal  curvature  is  zero,  the  corresponding  axis  will  be  at  infinity. 
Obviously  rotations  of  these  three  types  span  ker  dgi\). 
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Fig.  3.  Situation  where  the  contact  point  has  only  one  infinitesimal  degree  of 
freedom  (cf.  Section  5.2):  front  view  and  lateral  view. 

5.2.  Multiply  constrained  motion 

Now  we  consider  the  configuration  manifold  C fl  K of  a motion  constrained 
by  surface-surface  contact  and  a line  congruence  JC  as  described  in  Sect.  4.4. 
The  contact  point  still  depends  smoothly  on  g G C fl  K.  But  does  it  have 
nonzero  velocity  for  all  infinitesimal  motions  different  from  the  rotations  about 
the  contact  normal?  Obviously,  that  depends  on  the  intersection  ker  dgijj  fl 
Tg(C  fl  K).  Because  ker  dgip  C TgC,  this  intersection  equals  ker  dgip  fl  TgK. 
Depending  on  its  dimension,  there  are  the  following  three  possibilities: 

• dim  = 1:  Only  the  infinitesimal  rotations  about  the  contact  normal  are  in 
ker  dgip.  The  infinitesimal  motion  of  the  contact  point  is  two-dimensional. 

• dim  = 2:  The  rank  of  ip  is  one,  and  the  contact  point  varies  infinitesimally 
only  in  one  direction. 

• dim  = 3:  rk-0  = 0 and  the  contact  point  does  not  move  infinitesimally. 

If  the  line  g(a)  belonging  to  the  current  position  g € CC\K  is  a hyperbolic  line 
of  1C,  then  the  rank  of  ip's  restriction  to  CC\K  can  be  determined  geometrically 
(see  Fig.  3): 

Proposition.  We  use  the  notation  of  the  propositions  in  Sect.  4.3  and  of  the 
lemma  in  Sect.  5.1.  If  there  are  lines  S\,S2  such  that  Si  intersects  both  b1 , b" , 
Si  is  incident  with  p+n/ Ki  (i  = 1,2),  and  si  C p+[n]  + [w2],  S2  C p+[n]  + [wi], 
then  vkip\C  OK  < 1,  otherwise  rktp\C  fl  K — 2. 

Proof:  The  lemma  in  Sect.  5.1  describes  ker9  ip  A and  the  first  proposition 
in  Sect.  4.2  does  the  same  for  TgK  =:  B.  To  compute  A fl  B we  apply  the 
duality  which  assigns  to  a linear  space  A C 1R6  the  linear  space  A*  of  those 
infinitesimal  motions  (a*,  a*)  which  fulfill  a • o*  + a • a*  = 0 for  all  (o,  o)  € A. 
Clearly,  dim  A*  = 6 — dim  .<4  and  ( A fl  B)*  is  the  linear  span  of  A* , B*. 

If  is  well  known  that  an  infinitesimal  rotation  is  in  A * if  and  only  if  its  axis 
intersects  all  the  axes  of  the  infinitesimal  rotations  of  A.  Thus,  the  rotation 
axes  in  A*  consist  of  two  pencils  with  the  same  vertices  as  those  of  A,  but 
with  orthogonal  planes;  and  the  rotation  axes  of  B*  are  just  the  lines  6',  b" . 
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Rank  < 1 means  dim(A*+.B*)  < 4,  or  that  the  span  of  b'  and  A*  contains 
b".  It  is  easily  seen  that  apart  from  even  more  special  cases,  this  happens  if 
the  axes  of  this  span  are  a hyperbolic  linear  congruence  with  axes  si,  S2-  □ 

§6.  Collision  Checking 

It  is  possible  that  a milling  tool  has  no  local  intersections  with  the  finished 
surface,  but  while  in  contact  with  the  surface  at  one  point,  it  cuts  into  another 
part.  Algorithms  which  test  for  this  type  of  intersection  of  two  bodies  in  space 
are  time  consuming,  and  therefore  we  want  to  circumvent  the  general  collision 
test  in  some  way.  In  some  cases  we  are  able  to  predict  the  total  absence  of 
collisions  based  only  on  the  curvatures  of  the  boundaries  of  the  two  bodies 
involved. 

We  say  that  a surface  $ is  millable  by  a body  E,  if  (a)  there  are  no  local 
intersections,  and  (b)  it  is  possible  for  E to  move,  within  previously  imposed 
constraints,  along  the  surface  $ such  that  it  touches  $ during  this  motion, 
but  never  actually  intersects  it. 

6.1.  3-axis  milling 

3-axis  milling  means  that  a milling  tool  E rotates  about  its  axis  a and  moves 
such  that  a remains  parallel  to  a fixed  line,  and  E always  touches  the  finished 
surface.  As  the  rotation  about  a is  not  important  for  geometric  considerations, 
we  disregard  it  completely  and  consider  a body  which  moves  in  a translational 
manner. 

We  assume  that  E as  well  as  the  workpiece  $ are  bounded  by  piecewise 
twice  continuously  differentiable  surfaces  (convex  edges  are  allowed).  Then 
the  so-called  general  offset  surface  (defined  below,  see  Fig.  4,  right)  of  $ with 
respect  to  E shows  in  its  singularities  and  self-intersections  the  singular  posi- 
tions of  the  motion  and  the  collisions  (cf.  [15,16]). 

Definition.  Choose  a reference  point  p attached  to  E.  Consider  the  set  of 
translations  t such  that  r(E)  touches  (disregarding  intersections).  Then 
the  set  of  all  points  r(p ) is  called  the  general  offset  surface  of  $ with  respect 
to  E. 

Theorem  7.  If  E is  strictly  convex  with  positive  principal  curvatures,  and  <3? 
is  connected,  then  the  general  offset  surface  of  4>  with  respect  to  E is  smooth. 
It  is  regular  in  all  points  which  correspond  to  translationally  regular  contact 
positions,  and  if  it  is  both  regular  and  free  of  self-intersections,  then  $ is 
globally  millable  by  E. 

Proof:  (cf.  [16,23])  For  all  points  p € $ there  is  a unique  position  g such  that 
g( E)  touches  $ in  p.  The  parametrization  of  the  general  offset  which  thus 
is  induced  by  the  parametrization  of  $ is  easily  seen  to  be  smooth.  There  is 
always  a point  po  such  that  E D $ consists  of  po  only.  If  E touches  in  p,  and 
thereby  cuts  into  another  part  of  4>,  then  let  E move  such  that  the  contact 
point  follows  a curve  which  joins  po  with  p.  At  the  first  time  that  the  set 
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Fig.  4.  Left:  Line  congruence  /C,  finished  workpiece  $ and  cutter  E.  Right: 
Workpiece  4>,  milling-tool  E and  general  offset  surface  I\ 

E fl  $ consists  of  more  than  one  point,  E and  $ touch  each  other  in  all  points 
of  E fl  which  leads  to  a multiple  point  of  the  general  offset.  □ 

Thus,  we  have  transformed  the  collision  problem  into  the  problem  of 
determining  singularities  and  self-intersections  of  certain  surfaces.  If  4>  has 
some  additional  properties  then  it  is  not  difficult  to  guarantee  total  absence 
of  unwanted  collisions  provided  that  no  local  collisions  occur: 

• if  $ is  convex  and  E,  $ are  oppositely  oriented  (i.e.,  the  body  bounded 
by  $ is  its  inside);  or 

• if  $ is  convex  and  E,  $ are  equally  oriented  (i.e.,  the  body  bounded  by 
$ is  its  outside  and  E is  inside);  or 

• if  $ is  star-shaped  and  bounds  its  inside;  or 

• if  $ is  star-shaped,  bounds  its  outside,  and  E fits  into  the  convex  core  of 
$ (the  convex  core  of  a star-shaped  set  M is  the  set  of  points  with  respect 
to  which  M is  star-shaped  — it  is  a convex  subset  of  M);  or 

• if  $ is  the  graph  surface  of  a function  over  a planar  domain  whose  bound- 
ary is  millable  by  the  ‘top  view’  contour  of  E. 

Proof:  The  proof  of  this  can  be  found  in  [16],  and  generalizations  actually 
unimportant  for  applications  are  studied  in  [23].  The  idea  of  the  proof,  which 
is  important  also  for  the  proof  of  Th.  8,  is  as  follows:  We  assume  that  there 
is  a ‘projection’  of  entire  space  onto  4>.  If  e.g.  $ is  convex,  just  choose  any 
interior  point  o,  and  to  project  a point  p,  intersect  the  ray  op  with  4>. 

Then  consider  the  following  mapping  whose  domain  is  $:  For  a point 
x € $,  translate  E such  that  it  touches  $ there,  and  project  a previously 
chosen  reference  point  of  E onto  4>.  This  mapping  is  shown  to  be  smooth 
and  orientation-preserving,  and  by  scaling  E with  a factor  A (1  > A > 0) 
is  deformed  into  the  identity  mapping.  Differential  topology  allows  now  to 
conclude  the  $ is  one-to-one  and  onto,  which  means  that  no  translate  of  E 
can  touch  $ in  more  than  one  point.  □ 

6.2.  5-axis  milling  constrained  by  a line  congruence 

Here  5-axis  milling  means  that  the  milling  tool  moves  such  that  it  touches  the 
finished  surface,  and  its  axis  is  always  contained  in  a line  congruence  1C. 
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Definition.  A regularity  domain  of  1C  is  an  open  domain  M in  ]R3  such  that 
K.  defines  a fibration  of  M (every  point  of  M is  contained  in  exactly  one  set 
l fl  M,  where  l is  a line  of  the  congruence).  A section  of  a regularity  domain 
M is  a surface  which  intersects  all  sets  l fl  M exactly  once.  The  regularity 
domain  M is  a tubular  neighbourhood  of  its  section  $ if  it  is  diffeomorphic 
to  $ x ]R. 

The  domain  M is  a tubular  neighbourhood  of  if  the  sets  l fl  M are  open  line 
segments,  the  point  p = $ fl  l fl  M is  an  interior  point  of  this  segment,  and 
the  initial  and  end  points  of  the  segment  depend  smoothly  on  p. 

We  say  that  a convex  body  E with  rotational  symmetry  is  admissible  for 
a line  congruence  1C  and  a connected  closed  surface  <f>  which  is  the  boundary 
of  a subset  of  ]R3  if  the  following  is  fulfilled: 

• There  is  a regularity  domain  M for  K,  which  is  a tubular  neighbourhood 
of  $. 

• E moves  such  that  it  is  entirely  contained  in  the  regularity  domain. 

• In  no  position  of  C fl  K the  contact  tangent  plane  is  parallel  to  the  axis 
of  E. 

• The  contact  point  has  two  infinitesimal  degrees  of  freedom  for  all  scaled 
versions  AE,  0 < A < 1 (cf.  Sect.  5.2). 

These  conditions  are  actually  easy  to  fulfill  in  practice  except  for  the  last  one, 
which  is  difficult  to  detect  in  advance.  In  the  special  case  of  three-axis  milling 
(/C  is  a bundle  of  parallel  lines),  this  requires  that  all  contact  points  on  E are 
elliptic  surface  points. 

Suppose  that  we  are  given  a surface  $ and  a milling-tool  E,  and  we  have 
chosen  a congruence  1C,  and  have  found,  for  all  contact  points  on  <3>,  a position 
g( E)  such  that  p(E)’s  axis  is  in  K.  Suppose  we  have  already  tested  for  local 
millability  and  the  admissibility  conditions  described  above.  Then  we  have 
the  following 

Theorem  8.  Under  the  circumstances  described  above,  the  cutter  does  not 
interfere  with  the  surface  E,  i.e.,  $ is  globally  millable  by  E. 

Note  that  in  many  cases  it  will  be  sufficient  to  check  the  admissibility  only  for 
the  cutter  head,  because  collisions  of  the  cutter  shaft  with  the  workpiece  will 
be  treated  by  different  methods  (see  the  introduction) 

Proof:  In  Section  5 we  have  established  that  the  contact  point  depends 
smoothly  on  E’s  position  g.  Let  a be  E’s  axis,  and  let  p be  a rotation  about 
g(a).  Clearly  p o g is  again  a contact  position  with  the  same  contact  point  as 
g.  If  we  choose  a reference  point  p on  E’s  axis  and  inside  E,  the  path  surface 
of  p also  depends  smoothly  on  g,  and  p o g(p)  = g(p). 

This  means  that  the  contact  point  depends  smoothly  on  the  position 
g(p)  of  the  reference  point.  The  last  admissibility  condition  ensures  that  also 
g(p)  depends  smoothly  on  the  contact  point.  Thus,  we  can  define  a smooth 
mapping  / : $ — > $ as  follows:  A contact  point  q e $ is  mapped  to  the 
corresponding  point  g(p),  which  is  subsequently  mapped  to  the  intersection 
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f(q)  of  g(a)  with  4>.  If  we  scale  E with  a factor  A (0  < A < 1),  we  get 
mappings  f\,  where  fo  is  the  identity  mapping  of  $ onto  itself,  and  f\ (q) 
depends  continously  on  A.  f\  is  never  singular  since  we  are  in  a regularity 
domain,  /i  = / is  orientation-preserving  because  fo  is,  and  the  number  of 
pre-images  of  a point  is  the  same  for  fo  = id  and  f\  = /.  This  shows  that 
/ is  one-to-one  and  onto.  If  <?(E)  touches  $ in  one  point  and  cuts  into  $ in 
another,  there  also  is  a position  g( E)  where  E touches  in  two  different  points 
(cf.  the  proof  of  Th.  7).  These  two  points  have  the  same  /-image  by  definition 
of  /,  which  contradicts  bijectivity  of  /.  □ 
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